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The chemistry of extragalactic carbon stars 
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ABSTRACT 

Prompted by the ongoing interest in Spitzer Infrared Spectrometer spectra of carbon 
stars in the Large Magellanic Cloud, we have investigated the circumstellar chem- 
istry of carbon stars in low-metallicity environments. Consistent with observations, 
our models show that acetylene is particularly abundant in the inner regions of low 
metallicity carbon-rich AGB stars - more abundant than carbon monoxide. As a 
consequence, larger hydrocarbons have higher abundances at the metallicities of the 
Magellanic Clouds than in stars with solar metallicity. We also find the oxygen and 
nitrogen chemistry is suppressed at lower metallicity, as expected. Finally, we calculate 
molecular line emission from carbon stars in the Large and Small Magellanic Cloud 
and find that several molecules should be readily detectable with the Atacama Large 
Millimeter Array at Full Science operations. 
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1 INTRODUCTION 

The advent of the Atacama Large Millimeter Array (ALMA) 
and other new (sub-)millimetre telescopes, with their un- 
precendented spatial resolution and sensitivity, will allow 
the observation of giant stars in other galaxies in similar de- 
tail to that achieved for Galactic objects. These advanced 
capabilities prompt investigation into the nature of these 
extragalactic stars, with one of the most interesting aspects 
being the study of the effect of sub-solar metallicities on 
circumstellar chemistry and dust composition. 

Recent studies in the infrared using the Spitzer Space 
Telescope and ground-based instruments have highlighted 
the deep molecular absorptions of, primarily, acetylene in 
the spectra of evolved carbon stars, in the Magellanic 



The Magellanic Clouds (MCs) are nearby dwarf 
galaxies with sub-solar metallicities. The Larg e Magel- 
lanic Cloud (LMC), at a distance of ~50kpc (|Schaeferl 
|200S| ) . has an average metallicity of around 50% solar 
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absorption features are in general deeper than those seen in 
Galactic stars, and this implies that there is a difference be- 
tween the chemistry of Magellanic circumstellar envelopes 
(CSEs) and Galactic carbon stars, which are comparatively 
well-studied. 
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2008), has an average metallic it y ol around oO/o solar 
i Dufour. Shields, fe Talbotl Il982l : (westerlundl Il997ft . The 
Small Ma gellanic Cloud (SMC) , at a slightly larger distance 
of 66 kpc (fSzewczvk et~al]|2009h . has an average m et allicity 
of 20% solar <|Dufour et al.lll982l : IWesterlundlll997l 1. The ef- 
fect of this low metallicity regime on dust composition and 
galactic dust budg ets has been studied observationally by 
many authors (e.g..lZiilstra et al.ll2006l; lLagadec et alll2007i : 
van Loon et al.ll2008l : iMatsuura et all [2009). However, the 
effects on circumstellar chemistry have not been, as yet, 
studied in any detail, and our aim here is to pioneer in the 
field with this work. 

The chemical modelling of Galactic AGB stars was 
first attempted 35 years ago. Initial attempts focused 
on simple physical models and chemistry appropriate 
for oxy gen-rich (n(O) > n(C)) circumstellar environ- 
ments dGoldreich fc ; ScovnTel 1 19761 : IScalo fc Slavskvl Il980l : 
lJura fc Morris! Il98ll ). Physical models have largely re- 
mained simple (with some notable exceptions, e.g., 
ICordiner fc Millarll2009l l whereas the chemical modelling has 
moved to focus on carbon-rich {n(C) > n(0)) circumstellar 
chemistry since it shows a wider variety of molecules (e.g., 
Irluggins fc Glassgold|[l982l ). Progress in chemical modelling 
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is driven in part by the desire to explain observed abun- 
dances of newly-detected molecules in the most accessible 
carbon star, IRC+10216, for instance, in the a ddition of 
long ca rbon-chain molecu les (I Millar et al.l I2000T ) or anion 
species jMillar et al.ll2007T l. In this case, advances in technol- 
ogy have driven us to investigate carbon-rich circumstellar 
chemistry in previously challenging locations. 

In this paper, we present the results of modelling the 
circumstellar chemistry around carbon-rich AGB stars at 
sub-solar metallicities. We focus on three fiducial models, 
at metallicities and with physical conditions appropriate for 
the Galaxy, the LMC and the SMC. We initially assume so- 
lar metallicity (Z=0.02) for Galactic carbon stars, and aver- 
age interstellar metallicities for LMC (Z=0.008) and SMC 
(Z=0.004) carbon stars. In §2 we describe our adopted phys- 
ical model and discuss the different physical and chemical 
considerations, and also discuss sources of uncertainties in 
our models in §2.6. In §3 we describe how we calculate the 
chemical evolution of the circumstellar envelope. We show 
our results in §4, preceding a discussion in §5 (including our 
calculations of molecular line emission which may be observ- 
able with ALMA) and finally, in §6, we draw our conclusions 
regarding the chemistry of low-metallicity carbon stars. 



2 CHEMICAL AND PHYSICAL 

CONSIDERATIONS AT LOW METALLICITY 

Nucleosynthetic products dredged (via convection) from the 
stellar interior are mixed to the surface of the star, where 
material is accelerated to the ter minal velocity of the stel- 
lar w ind within a radius of 20 R* (Ke adv. Hall. fc Ridgwavl 
Il988f ) and passed into the CSE. The gas, which is mainly 
molecular hydrogen, is well-mixed with dust grains. We as- 
sume a spherical geometry where the gas has a 1 / R 2 density 
distribution, and a temperature profile which follows 

T(r) = max[150(E/Ro)" ' 79 ; 10] K. (7? ^ R = 5xl0 15 cm) 

(e.g.. iMillar et all l2000l : iMillar fc Herbstl Il994h . The CSE 
is irradiated by the interstellar radiation field (ISRF), 
but not by UV photons from the star itself, which are 
quenched in the stellar atmosphere. Extinct ion in the CSE 
is calc ulated according to the approach of ljura fc Morris! 
(|l98ll ). assuming interstella r-type grains, and we treat CO 
self-sh ielding according to M amon. Glassgold. fc Hug gins 
(1988). We model the carbon-rich chemistry in the circum- 
stellar region between ~200-100 000 R* (0.005-3 x 10 18 cm). 

In the subsequent sections, we consider the chemical 
and further physical ingredients of our model. Observation- 
ally, very little is known about the chemistry of extragalac- 
tic carbon stars, and so we mainly discuss the chemistry of 
Galactic carbo n stars, using the nucleosynthesis models of 
iKarakasl (|201Qh to adjust for lower metallicities (see H2.ip . 
Physical parameters of LMC and SMC carbon stars are more 
well constrained observationally, and we summarise those 
aspects in 32T2TJ231 

2.1 Initial chemical abundances 

As circumstellar material cools during the expansion of the 
stellar wind, it is energetically favourable for atomic ele- 
ments to form molecules in the photosphere of the star 



(R ~5 5R*), where the chemistry is in thermal equilib- 
rium (TE) due to the high temperatures and densities. For 
molecules of high stability, that is to say, those which con- 
tain strong internal bonds, the abundances amounted in 
this region of the star are carried through into the non-TE 
outer envelope (R ^ 100 R*) without significant change. In 
carbon-rich environm ents such molecules inc lude CO, N2, 
HCN and C 2 H 2 (e.g.. lTsuii|[l973l ; lMillarll2008h . 

A further group of molecules are also observed at high 
abundances at the inner edge of the outer envelope in Galac- 
tic stars. This subset is generally formed at low abundance 
in the TE photosphere of the star, but attains a higher abun- 
dance in the inner envelope where dust is formed and then 
accelerated (5 is 7? ?S 100 R*) due to non-equilibrium pro- 
cesse s. They comprise CS, SiO, SiS, NH3, SiEU and CH4 
(e.g., ICherchnefj[200r3 ). These two subsets of molecules are 
often called "parent species". 

To obtain accurate initial abundances for these ten par- 
ent species when modelling the chemistry in the outer en- 
velope (R i; 100 R*) at different metallicities, one must 
employ different techniques for the two subsets of par- 
ents. For the high stability parent species we can use 
a TE model to calculate the relevant data. Thi s ap- 
proach has been used previous l y to good effect (e.g., iTsuiil 
1 19731 ; ISharp fc Huebneri Il99d ; IMarkwickl |2000| ). For the 
remaining parent species other physical factors must be 
taken into considera tion, such as depletion through dust 
formation (SiS, Si O; iBieging fc Nguven-Quang-R icu 1989; 
Schoicr et al. 2006) , pulsati on-driven shocks (SiS, SiO, CS; 
Willacy fc Cherchnefllll998r ) and gas-grain interaction (e.g., 
hydrogcnation leading to NH3, S1H4 and CH4). Modelling 
these complexities in detail, which indeed themselves have 
been the subject of much investigation, is too advanced for 
this initial study and thus we must use other arguments, 
discussed below. 

2.1.1 TE molecules (CO, N 2 , C 2 H 2 , HCN) 

The TE model used is similar to that detailed by 
ISharp fc Huebneri (|l990T ) and lMarkwickl (doOOjE 

which work 

by minimising the Gibbs free energy of the thermodynamic 
system. We adopt a temperature of 2 250 K and a pres- 
sure of 1.033xl0 -3 atm, a ppropriate for the pho t osphere of 
a car bon-rich AGB star ((ivezic fc Elitzurl 1 19961 ; IMarkwickl 
2000). Results of the TE calculations for the range of tem- 
peratures 2 500-2 000 K are displayed in Fig. [TJ those species 
for which TE is appropriate are largely invariant across the 
temperature range. Elemental abundances for AGB stars at 
different me tallicities are ta ken from the nucleosynthesis cal- 
culations of Karakas] (|20 lOh for a 3Mq star. This mass was 
chosen because it was the lowest-mass model for which the 
star became carbon-rich at the three metallicites Z=0.02 
(Galaxy), Z=0.008 (LMC) and Z=0.004 (SMC). There are 
incongruities with these models, in that it is clear from ob- 
servational evidence that carbon stars in the Milky Way can 
form at masses as low as ~1.5M0 i Wallerstein fc Knappl 
Il998l ). The se modelling is sues are acknowledged and being 
addressed (|Karakadl201ll ), and we continue to use the re- 
sults since they are a consistent set of readily-available abun- 

1 Available on request to the authors. 
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Table 1. Initial model fractional abundances (w.r.t. H2) 



opGCIGS 


Observed Lralactic inner 
envelope abundance 


± Hi abundance 
(Galactic) 


Adopted initial 
abundance (Galactic) 


Adopted initial 
abundance (LMC) 


. . . 

Adopted initial 

abundance (SMC) 


CO 


a 6xl(T 4 


4.4xl0~ 4 


4.4X1CT 4 


1.7X1CT 4 


8.9X1CT 5 


N 2 


b 


5.4X1CT 5 


5.4xl0- 5 


1.2xl0~ 5 


3.9X1CT 6 


C2H2 


c 8xl0" 5 


3.3xl0~ 5 


3.3xl0~ 5 


2.8xl0" 4 


3.1 xl0~ 4 


HCN 


d 2xl0~ 5 


2.4xl0" 5 


2.4xl0~ 5 


3.3xl0~ 5 


2.0xl0~ 5 


SiS 


e 3xl0~ 6 


3.7xl0- 6 


6.9xl0~ 6 


2.1xl0~ 6 


l.OxlO- 6 


CS 


c 4xl0~ 6 


4.7xl0~ 6 


2.4X1CT 6 


1.6xl0~ 6 


8.0xl0~ 7 


NH 3 


^'2xlO~ B 


3.3X10- 11 


2.0xl0~ 6 


8.0xl0~ 7 


4.0xl0" 7 


CH 4 


c 2xl0" 6 


6.1xl0~ 9 


2.0xl0~ 6 


8.0xl0~ 7 


4.0xl0~ 7 


SiO 


9 lxl0~ 6 


2.0xl0~ 8 


1.7xl0~ 6 


6.8xl(T 7 


3.6xl0~ 7 


SiH 4 


c 2xl0~ 7 


8.8x10" 15 


2.0xl0~ 7 


8.0xl0~ 8 


4.0x10-* 



"iKwan fc Linkel il982l) ''Not observed due to lack of permanent dipole. 2x10 
ICordiner fc Millarl 120091) . similar to the elemental abundance ^Keadv 



iMonnier et al 



.. Ridgwavl 

J2007a|l: ICernicharo et all lll999h HSchoier et al.l ll2007bh: iDecin et al] , 
l|2000l ); lKeadv fc Ridgwavl il993l 'i I ^Schoier et al.l ||2006| 1 ; iKeadv fc Ridgwavl dl993h 



u sually assumed in mode ls (e.g. . 
^993); ICernicharo et al I lll999ll 
201ol ) - iHasegawa et al] (2006li : 



dances. Average stellar yields were used as inputs for the TE 
model and the resulting fractional abundances with respect 
to H 2 can be found in Table[T]for CO, N 2 , C 2 H 2 and HCN. 
This method (using the results of nucleosynthesis calcula- 
tions), is preferable to using the elemental abundances of 
the Magellanic Clouds in general, since much of the circum- 
stellar chemistry depends on the a mount of carbon gener - 
ated by the star in the AGB phase dMatsuura et al.l l2008). 
In the next section, we describe our method of determining 
input abundances for those molecules for which TE does not 
apply. 



2.1.2 Non-TE molecules (CS, SiO, SiS, NH 3 , SiH 4 , CH 4 ) 

The observed distribution of SiO peaks within a few 
stellar radii in the archetypal carbon star , IRC+10216 
l|Schoier et al l l200rj ; IKeadv fc Ridgwavl 1 19931 ). SiO is seen 
at high fractional abundances, n(SiO)/n(H 2 )~10 -6 , whereas 
TE models underpredict the abundance by around a factor 
of 30. Its formation is thought to be due to circumstellar 
shocks, and hence its abundance is dependent on that of 
atomic oxygen, since in the shocked material the reaction: 



Si +OH - 
dominates 



SiO +H 



(IWillacv fc Cherchnefll 



1980). 



(1) 

1 19981 ; 

As such, 
of SiO to 



Irlartquist. Dalgarno. fc Oppenheimei 
we choose to set the fractional abundance 
10~ 3 n(O), in line with that observed in IRC+10216. 

SiS is the main repository for both silicon and sulphur 
in the gas phase, and forms at a high abundance in TE 
models. Observations indicate that it is rapidly depleted as 
gas flows through the intermediate envelope, presumably 
due to the adsorption of the molecule onto grain surface s 
(|Bieging fc Nguven-Quang-Rieul Il989l ; I Boyle et all Il994h . 
Thus we assume that the initial abundance of SiS is equal 
to the elemental abundance of sulphur minus the fraction of 
CS (discussed below). For IRC+10216 this premise agrees 
well with observations jBieging fc Tafallalll993l : II ucas et al.l 
Il995l ; ISchoier et al.ll2006l ; IDecin et al.ll2010t ). 

CS is slightl y overproduced in TE m odels compared 
to observations (IKeadv fc Ridgwavl [l993) and work by 
IWillacv fc Cherchnefll (| 19981 ) has shown that CS is destroyed 



by pulsation-driven shocks by a considerable amount, the ex- 
tent of which depends on the assumed shock speed. It then 
slowly reforms in the outflowing wind. We assume slow shock 
sp eeds (lOkms ) and the sam e shock-destruction efficiency 
as I Willacy fc Cherchnefll l|l998l ). and adopt an abundance of 
CS a factor of 2 less than its TE abundance. 

Emission from the hydrogenated species NH3, CH4 
and SiEU originates in the dus t-formation zone and in - 
ner envelope around IRC+10216 (|Keadv fc Ridgwavlll993l ). 
Although surprising because of the high temperature of 
the dust in these regions, these molecules are most likely 
formed via hydrogen-addition reactions on the surface of 
dust grains. The very short time scale before ejection from 
the grain surface would indicate that hydrogen is the dom- 
inant reaction partner for adsorbed atoms. Typically these 
molecules a re of low abundance in both TE models (e.g., 
Tsuiil 1 19641 ). and non- TE models (|Willacv fc Cherchnefll 
19981 . and others), adding further weight to grain-surface 
hydrogenation theories. Observations of warm NH3 in 
IRC+10216 have been used to derive fractional abundances 
of 2-2.0xl0~ 6 llKeadv fc Ridgwavl 1 19931 ; IMonnier et all 
l2000l ; lHasegawa et al.ll2006l ). We assume that the abundance 
of NH3 at other metallicites will scale with the dust surface 
area. For the LMC and SMC this quantity, and the associ- 
ated gas-to-dust ratio and grain-size distribution, are highly 
uncertain, although work is currently ongoing to determine 
more definite values. Assuming that the grain-size distribu- 
tion in LMC and SMC carbon stars is similar to that in 
Galactic carbon stars, and that the gas-to-dust ratio scales 
with metallicity, we will adopt the upper observed fractional 
abundance of NH3 for Galactic carbon stars, and scale ac- 
cording to metallicity for the LMC and SMC carbon stars. 
We apply a similar reasoning for scaling the initial fractional 
abundances of CH4 and SiH.4 with metallicity. 

Adopting the initial abundances detailed in Table [1] 
means that only small amounts of nitrogen, oxygen and sul- 
phur are available for incorporation into other molecules and 
dust. However, in the Galactic case, some 5% of the ele- 
mental carbon remains for incorporation into molecules or 
carbonaceous dust. At LMC and SMC metallicities, the per- 
centage of free carbon rises to 28% and 32%, respectively, 
which is in reasonable agreement with iFerrarotti fc Gaill 
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Figure 1. Fractional abundances (log scale) at thermal equi- 
librium for MW, LMC and SMC metallicities, at a pressure of 
P~10 — 3 atm. We utilise values at 2 250 K in the subsequent mod- 
elling. 



(2006) for an evolved star. These authors also found that 
at most 9% of elemental silicon condensed into dust. Our 
initial abundances mean that 52-60% of the elemental sili- 
con is available for dust or dust seed formation. 

In Table [2] we have compiled a comparison of C/O ra- 
tios and the carbon excesses (log(ec — eo)) for the three 
metallicity regimes. C/O ratio s increas e dras tically in the 
nucleosynthesis calculations of Karakas (2010), largely due 
to the decrease in elemental oxygen. The carbon excess is 
large r at low metallicity. an effect which has been observed 
(e.g., IWahlin et ai1l2006n . 



Table 2. C/O ratios, elemental carbon and oxygen abundances 
( e C; e o)i an d carbon excesses in the three metallicity environ- 
ments. 



Environment 






C/O 


C-O 




Milky Way 


9.05 


8.94 


1.3 


2.5x10" 


-4 


LMC 


9.33 


8.53 


6.3 


1.8x10" 


-3 


SMC 


9.33 


8.25 


12.0 


2.0x10" 


-3 



2.2 Mass-loss rate 

Mass-loss rates for AGB stars are generally measured using 
observations of the CO envelope and then assuming a CO- 
to-Hb scaling ratio (the so-called 'X-factor', Xco)- Alterna- 
tively, the near- and mid-infrared spectral energy distribu- 
tion can be fitted with a radiative-transfer model to obtain 
a dust mass-loss rate, and then by assuming a gas-to-dust 
ratio, the total (gas+dust) mass-loss rate can be calculated. 
The former method has been routinely used for Galactic 
stars, but extragalactic stars are in general too faint. The 
availability of high-sensitivity data from, for example, the 
Infrared Space Observatory or the Spitzer Space Telescope, 
means that the second method is more practicable for extra- 
galactic sources, and a large nu mber of authors have ta ken 
this approach - see for ex ample Ivan Loon et al.l l|l999rj) or 
ICroeneweeen et al.l (|2009l ) and references therein for mass- 
loss derivations for stars in the LMC and SMC. However, 
these two met hods probe different parts of the circumstellar 
envelope (e.g.. iKemper et al.l [2003). and often give differing 
values of total mass-loss rate. 

Mass-loss rates for AGB stars in the Galaxy range 
from 10 -9 to 10~ 4 Mp ) yr~ 1 , with the median bein g 
3xlO- 7 M yr- 1 l|Olofssonll2008bl ; ISchoier fc OlofssorfeOOll ). 
A typical mass-loss rate for a carbon star with a rich cir- 



cums tellar chemistry suc h as IRC+10216 is ~10 M^yr 
(e.g.. IWoods et al.l l2003i : iMen'shchikov et al.ll200ll ). There 
seems to be no differentiation in mass-loss rate based on 
chemistry (M-, C- and S-stars) in the Galactic sample of 
lOlofssonf (|2008bl ) (see Fig. 2). There is a weak dependence 
on metallicity for O-rich stars compared to C-rich stars, with 
the difference ap pearing to be due to the driving mechanism 
of the wind (e.g. Hofncr 2007), and the opacity of the dust 
grains within it (Woitkc 200<I). 

In the halo of our galaxy, where the metallicity is lower 
than that in the plane, a sample of 16 carbon-rich AGB 
stars have mass-l oss rates of ^ 4xl0~ 6 Mpyr -1 , to within 
a factor of three |Mauron 20 08|). Similar, al t hough smaller, 
rates are estimated by Groenewegen et all (|l997l ) for two 
of these halo stars. In the metal-rich - but still sub-solar 
- Sagittariu s dwarf spheroidal ga laxy, six carbon stars are 
detected by lLagadec et al.l (|2010l ) , who estimate mass- loss 
rates on the order of 10 -6 M©yr _1 . 

In the LMC and SMC, carbon st ars have simila r mass- 
loss rates to those in the Galaxy. iTanabe et all (|l997l ) 
found that dusty carbon-ric h LMC stars have mass-loss 
rates less than 10~ 5 Mmvr" 1 ; iLeisenring. Kemper, fc Sloanl 
(2008) found that the brightest Magellanic carbon stars 
which made up their (biased) sample have mass-loss rates 
in a narrow range around lO _6 M0yr _1 ; Ivan Loon et alj 
(|2003h found a particularly high mass-loss-rate carbon 
star, LI-LMC 1813, which has a mass-loss rate of nearly 
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4x10 M0yr~ . Since this star is in a cluster, the birth 
mass and metallicity could be determined, allowing for 
an estimate of the dust-to-gas ratio (see i|2.3p . Vari- 
ous carbon stars have been observed in the SMC (e.g., 
van Loon. Ziilstra. fc Groenewegenl Il999al ; iMatsuura et al.l 
20051 ; I van Loon et al .1120081 ) and rates are found to be similar 
to those in the LMC (see also Ivan Loonll2000l , l2006l ). 

Since in this work we are only considering carbon stars, 
for our fiducial models in all three metallicity regimes we 
will assume M=3xl0~ 5 M Q yr _1 . This is toward the high 
end of the observed range of mass-loss rates, but typical 
of stars where the strong infrared C2H2 features are seen 
most clearly, and of the well-known Galactic carbon star 
IRC+10216. 



2.3 Gas-to-dust ratio 

In the early 1990s lHabing. Tignon. fc Tielens! (|l994h postu- 
lated that the dust-to-ga s ratio wou l d dep end on metallicity, 
and this was refined by Ivan LoonI ( 2000) who determined 
observationally that it had an approximately linear depen- 
dence. Dust production at low metallicity is limited by the 
availability of heavy metals (e.g., tit anium) to form the co n- 
densation seeds for dust formation l|van Loon et al.l l2008). 

There have been several attempts to determine gas-to- 
dust ratios in the diffuse ISM of the MCs from Hi and in- 
frared maps. At higher densities this becomes almost im- 
possible. Moreover, it is not clear how the gas-to-dust ra- 
tio in the diffuse ISM relates to that in circumstellar en- 
vironments. In the ISM of the LMC, the dust-to-gas ratio 
is approximately a quarter that of the solar value, and this 
has been established fo r several decades (e.g.. Ivan Gendered 
Il97d ; iKoornneej Il982l ; IClavton fc Martini 1 19851 ) and con- 
firmed more recently with the Herschel and Spitzer tele- 
scopes (e.g., iMeixner et ail l201fj| ; iGordon et atl l2003h . In 



that of solar (Ivan den Berghl 


1968; Ivan 


Genderen 


197C 


Leaueux et al.l 1982; 


Bouchet et al.l 


1 19851; Hot et all 


2004; 


Gordon et al.ll2009l). Polvcvclic Aromatic Hydrocarbons be- 



have differently, but a re also depleted at low metallicity 
jSandstrom et al.ll20l3 ). 

In the circumstellar envelopes of AGB stars, the dust- 
to-gas ratio is harder to determine, especially without a good 
measure of the gas mass. This ratio also depends on chem- 
istry: i n oxygen-rich envelo pes the dependence is roughly 
linear l|Marshall et al.ll2004l ). In carbon-rich envelopes there 
is some decline wi th metallicity, b ut possibly slightly shal- 
lower than linear (|van LoonI |2000J). Estimates deri ved from 
molecular band strengths (e.g., Sloan et all 120061 ) can be 
misleading (|van Loon et al.ll2008l ) . In lieu of firm determina- 
tions, authors have assumed various values: Ivan Loon et al.l 
|2003l ) assume a ga s-to-dust value of 300-50 for the carbon 
star LI-LMC 1813; iLeisenring et all |200Sl ) assume 100 for 
the Galaxy, 200 for the LMC and 500 for the SMC. We adopt 
100, 200 and 500, respectively, in line with the metallicity, 
for the fiducial models. 



2.4 Wind expansion velocity 

The expansion velocity of a circumstellar envelope, that is, 
the terminal velocity once the wind has undergone accelera- 
tion in and close to the dust formation zone, depends on the 



dust-to-gas ratio, which in turn depends on the metallicity. 
There is also a lesser reliance upon the luminosity of the star, 
such that Wexp oc i]/ 1 / 2 ^ 1 / 4 ^ where ^ is the dust-to-gas ratio, 
and L the luminosity dHabing et al.lll994l ; lElitzur fc Iverici 
l200ll ; iMarshall et alj|2004l ). Thus expansion velocity is ex- 
pected to be lower in the Magellanic Clouds, for a given 
luminosity. 

There is also a link between expansion velocity and the 
mass-loss mechanism, since pulsation-driven winds (which 
are associated with low mass-loss rates) show low expansion 
velocities. Stars experiencing a superwind, where the mass- 
loss driver is radiation pressure on dust grains, general ly 
show larger expansion velocities (e.g.. IWinters et al.l l2000). 

In general, expansion velocities in the LMC have been 
measured to be fairly small compared to the Galaxy. For 
all chemistries, Galactic expansion veloci ties cover a wide 
range, from low (e.g., 1.5-2 2.5 kms~ 1 : lOlofssonl l2008bl ) 
to high (e.g., 4.3-35.4 km s -1 ; [Loup et al Il993l )and even 
higher vel ocities are expected theoretically (e.g., up to 
eOkms -1 ; iMattsson. Wahlin. fc Hofnerl l2010l ). For Galac- 
tic carbon stars undergoing a superwind mass-loss, typ- 
ical expansion ve locities range from 13-22 km s" 1 ( e.g. , 
IWoods et aT]|2003l ). In the LMC. Ivan Loon et all |2003l ) de- 
rived «ex P =9.5kms _1 from modelling the carbon star LI- 
LMC 1813. In a sa mple of oxygen- r ich ci rcumstellar en- 
velopes in the LMC, [Marshall et alj (|2004h measured the 
expansion velocity of a number of OH masers and ob- 
tained results mostly in the region 8-17kms _1 . Circumstel- 
lar OH masers are thought to trace the terminal velocity 
of the wind, rather than H2O masers, for example, which 
trace the acceleration zone. IWood. Habing. fc McGregor! 
(|l99Sl ) suggested that stars with higher expans ion veloc- 
ities h ave higher metalliciti es. More recently, Ivan LoonI 
(2000); M arshall et alj |2004j, and references therein) and 
ISchoier et al.l ( 2007al ) have suggested that expansion veloci- 
ties increase with evolution on the AGB. Expansion veloci- 
ties for stars in the SMC have not been determined. It should 
be noted that gas velocities and dust velocities are not nec- 
essarily well coupled, and d rift velocitie s can be substantial 
at low mass-loss rates (e.g., Jones 2001). 

One can also look to the Galactic Halo as a low- 
metallicity environment. Expans ion velocities have b e en de - 
termined for several Halo stars: iGroenewegen et al.l (Il997l ) 
measured an expansion velocity of 3.2kms _1 from a 12 CO 
J=2-l line, which is exc eptionally low for a ca rbon star. This 
value was confirmed bv lLagadec et all l|2010l ). who observed 
the 12 CO J=3-2 line. The y also calculated a l ow (dust) 
mass- loss rate for this object. ILagadec et al.l (|2010l ) observed 
two other carbon stars in the Halo, and determined expan- 
sion velocities of 6.5 and 8.5 km s -1 . Three other carbon-rich 
objects in their sample were thought to be in the Galactic 
Halo, but instead reside within the Galaxy's more metal- 
rich thick disc. Expansion velocities for these stars were 
measured as 11.5-16.5 km s _1 . Thus it seems clear that for 
carbon stars, at least, expansion velocity decreases with de- 
creasing metallicity. This is consistent wi th theoretical mod- 
els of AGB outflows |Wachter et al.ll2008l ) , where low expan- 
sion velocities are due to the formation of less dust, which in 
turn implies a less efficient acceleration of the circumstellar 
material. 

For the purposes of our fiducial models, we have 
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Table 3. Physical parameters of the fiducial models 



Parameter 


MW 


LMC 


SMC 


M (Moyr- 1 ) 


3xl0" 5 


3xl0" 5 


3xl0- 5 


1/* 


100 


200 


500 


fcxp (kms" 1 ) 


20 


10 


5 


G (Go) 


1 


2 


4 



M is the mass-loss rate, \P the dust-to-gas ratio, i) e xp the envelope 
expansion velocity, and G the interstellar UV field. 

adopted u e xp=20, 10 and 5kms -1 for the Galaxy, LMC and 
SMC respectively. 



2.5 UV radiation field 



iBernard et alj (|2008l ) were able to establish that the in- 
terstellar radiation field in the LMC is in general stronger 
than that in the Milky Way. They find that it varies from 
0.8 Go in diffuse regions to 3.5 Go in molecular regions, with 
an average value of ^2 Go. We use this value in our mod- 
elling. Such a galactic variation in field strength is also seen 
in the Galaxy, where the strength at the inner molecular 
ring c an be up to five tim es that in the solar neighbour- 
hood dPaladini et all 120071 ). For the SMC, recent work by 
ISandstrom et al.l (|201dl . Fig. 10) has shown that the ISRF 
has an average value of ~4Go, although it can be as high 
as 10-30 Go in regions of star formation. 

2.6 Caveats 

Although we are presenting here three fiducial, and ex- 
ploratory, models, not intended to represent any particu- 
lar source, it is worth noting sources of uncertainties in our 
models, and how they may affect chemistry and predicted 
line intensities. Many of the physical parameters above have 
mainly an effect on the radius at which UV photons will 
dominate the chemistry of the envelope. Decreasing the 
mass-loss rate, increasing the expansion velocity of the wind 
or ISRF strength, or decreasing the dust-to-gas ratio (i.e., 
reducing the dust content) will have the effect of allowing 
greater penetration for UV photons, meaning that the en- 
velope chemistry transitions to a photochemistry (rather 
than an ion-molecule or neutral-neutral chemistry) closer to 
the star. In practice, this means that parent molecules are 
destroyed more effectively, daughter species are produced 
and then destroyed in a narrower shell, and the distribu- 
tion of species becomes more condensed. Since abundant 
parent species are found in the more dense parts of the en- 
velope, where column densities are high, any changes to this 
threshold radius will have a minimal effect on the emission 
intensity from parents. 




16.0 16.5 17.0 17.5 18.0 18.5 
log(R) [cm] 

Figure 2. Fractional abundances of TE parent species. The solid 
line indicates the Galactic model, the dotted line is the LMC 
model, and the dashed line is the the SMC model. 




16.0 16.5 17.0 17.5 18.0 18.5 
log(R) [cm] 

Figure 3. Fractional abundances of the polyyne species. Line 
characteristics are as in Fig. [2] 



481 chemical species, including atoms, neutral molecules, 
cations and anions, linked by 6174 reactions, which is consid- 
erably larger than the most recent release of UDfA, Rate06. 
Species range in size from single atoms (H, He, C, Mg, N, 
O, S, Si) to long hydrocarbon chains (e.g., C23, C23H2). 

We follow a parcel of gas as it passes from the inner 
edge of the circumstellar envelope outwards at the terminal 
wind velocity. Since we assume a constant mass-loss process, 
this procedure is time independent, and we obtain a snap- 
shot of the chemical structure of the circumstellar envelope 
during the AGB phase. Chemical rate equations are s olved 
using the Gear method for stiff differential equations |Gearl 
1 197 ll ). resulting in fractional abundances for all species at 
each radial gridpoint. 



3 CHEMICAL MODELLING 



Our chemical network is based on that o f lCordiner fc Millar! 
(2009), augmented with silicon chemistry fr om the UMIST 
Data base for Astrochemistry 2006 (UDfA; IWoodall et al.l 
2007). Reactions with Si-bearing species are included 
only for those species a lready present in the network of 
ICordiner fc Millar] (|2009l ). The resulting network contains 



4 RESULTS 

Radial abundance profiles for various species are plotted in 
Figs. [2H6] and we discuss different families of molecules as 
follows. 
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16.0 16.5 17.0 17.5 
log(R) [cm] 
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log(R) [cm] 



18.0 18.5 



Figure 5. Fractional abundances of the anion species. Line characteristics are as in Fig. [2] 




16.0 16.5 17.0 17.5 18.0 18.5 
log(R) [cm] 

Figure 4. Fractional abundances of HNC and cyanopolyyne 
species. Line characteristics are as in Fig. [2] 
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Figure 6. Fractional abundances of the silicon-bearing species. 
Line characteristics are as in Fig. [2] 



4.1 Carbon monoxide, polyynes and 
cyanopolyynes 

Carbon monoxide (CO) in Galactic carbon stars is the most 
abundant molecule after H2, and the most readily observ- 



able in the submillimetre regime, given its high abundance 
and significant dipole moment. However, it is clear from 
the calculations that at lower metallicity (even at just half 
solar metallicity in the LMC), acetylene (C2H2) becomes 
the dominant carbon-bearing molecule for 7? <10 16 cm 
(Fig. [2]). In terms of column density through the envelope, 
N(C 2 H 2 )>N(CO) for LMC metallicity and is three times 
greater at SMC metallicity. The CO abundance is lower in 
LMC and SMC carbon stars, and this reduces the efficiency 
of self-shielding, meaning that the extent of the CO envelope 
is smaller in comparison. 

The predominance of acetylene at lower metallicity 
means that the hydrocarbon chemistry is also more de- 
veloped in such environments (Fig. Acetylene and the 
ethynyl (C2H) radical are the basis of much of the carbon- 
chain growt h that occurs in c arbon-rich circumstellar en- 
vironments (|Millar et alj|200(ih . and so larger species such 
as triacetylene (C6H2) experience a large boost in pro- 
duction. In the LMC and SMC models, N(C 6 H 2 ) = 3.5- 
3.9x 10 17 cm -2 , whereas in the Galactic model the column 
density is 30 times smaller. The peak abundance of C 2 3H 2 , 
the largest polyyne in the model, is comparable for the two 
lower metallicity regimes, but ~100 times smaller for the 
Galactic model. This increase in large hydrocarbons also 
applies to the cyanopolyyne chains (HC3N, HC5N, etc.) de- 
spite the reduction in elemental nitrogen abundance (Fig. 3]). 
This would indicate that as metallicity decreases, nitrogen 
is preferentially sequestered in HCN and the cyanopolyynes 
rather than in other nitrogen-bearing species. Tabled which 
gives an overview of the nitrogen repositories in the model 
by comparing column densities through the envelope, shows 
that this is the case. The abundance of HdN (odd x ^ 3) 
species increases from Galactic metallicity to LMC metallic- 
ity but decreases at SMC metallicities, following the changes 
in initial fractional abundance of HCN (Table [TJ. 

The distribution of the cyanopo l yyne chains is very 
different to that seen in iMillar et all J2000, Fig. 6); how- 
ever it matches very well with the smoot h density distri- 
bution model of ICordiner fc Millar] l|2009h . The difference 
arises from the inclusion of an increased photodissociation 
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Table 4. Main nitrogen repositories at differing metallicities. 



Species 


MW 


LMC 


SMC 


N 2 


71% 


37% 


25% 


HCN 


14% 


42% 


58% 


N 


12% 


12% 


12% 


HC S (3...n)N 


1% 


6% 


3% 


NH 3 


1% 


1% 


1% 


CN 


<1% 


<1% 


<1% 



rate for HC3N, resulting in a substantial reduction of the 
HC3N/HC5N ratio. 

4.2 Anions 

Several anions have now been detected in Galactic carbon 
star I RC+10216: C 6 H" was identified by iMcCarthv et all 
(2006), and confirmed bv lCernicharo et al.l (2007), who also 
re ported the initial dete ction of C4H~. CsH - was de tected 
bv lRemiian et all (120071 ) and iKawaguchi et al.l (|2007l '). The 
anions of the cyanopolyyne radical s C3 N and C5N were 
discov ered by iThaddeus et all J2008l'l and Cernicharo et all 
(|2008t) . respectively. Most recentlv lAgundez et al.l (|20ld ') de- 
tected the smallest molecular anion to date, CN - . A sum- 
mary of observed column densities and estimates of neutral- 
to-anion ratios is sho wn in Table [5] The anion chemistry is 
discussed in detail bv lCordiner fc Millar! l|2009fl . 

In general, the radiative electron attachment 
rate of these carbon chains increases with length 
l|Herbst fc O samura 120081 ') . however, the abundance of 
the neutrals from which the anions are created peaks 
with CeH, and hence the most abundant anions in these 
models are CsH - and C7H - (Fig. [5]). We can compare 
the observed column densities and ratios for these anions 
with those calculated in our Galactic model (also shown 
in Table [5]). Despite not specifically modelling IRC+10216, 
the agreement between model and observation is reasonably 
good for the smaller hydrocarbons. For column densities 
of the larger hydrocarbon anions, C4H - and larger, the 
agreement is not as close, with the model overproducing by 
factors of 100-1000. The abundances of these hydrocarbon 
anions are particularly sensitive to t he initial abundance o f 
acetylene, as has been discussed by iRemiian et all (2007). 
They show that a reduction in the initial abundance of 
acetylene by a factor of 5-10 brings about a much better 
overall agreement. The C2H - anion has not yet been 
detected, and we predict a column density of ~10 9 cm -2 , 
several orders of magnitude below that of the detected 
anions. Thus, unless we have underestimated the radiative 
electron attac hment rate of C2H, C 2H~ is not likely to 
be detectable. lAgundez etHI (|2010h have determined an 
upper limit for this species in IRC+10216, of <0.0014% of 
the C2H column (Table [5]). Similarly, C10H - has not been 
detected, but has a similar ion-neutral ratio to CsH - and 
CgH~, and a column density which is about a quarter of 
that of CsH - , according to the model. Also abundant in 
the envelope are negatively-charged carbon chains, C7/ 9, 
which have column densities of 8x 10 13 -1 x 10 15 cm -2 . 

The distribution of hydrocarbon anions (Fig. [S| is sim- 
ilar for all members of the family from C4H - to CsH - . 
Similar distributions are also found for Ci,3,s,7N - . At the 



inner peak of its distribution, at log(7?) = 15.85, CN - 
is formed via electron a ttachment to M g NC, even in the 
lowest metallicity case. lAgundez et"al] (|2010h find that 
HCN + H - — y CN - + H 2 provides a minor contribution 
in this region, but we find this reaction ineffective. The 
much broader, outer peak of the CN - distribution arises 



due to the reaction N + C 4 



CN - + C 3 ...6- In gen- 



eral, the larger members of the family are all formed 
by electron attachment to corresponding neutrals. The 
distribution of CN - fits well with the observations of 
lAgundez et"al] (|2010h . and the column density calculated 
in the model (3.0xl0 12 cm -2 ) matches the observed deter- 
mination (8. Ox 10 12 cm -2 ) very closely. Similarly, the col- 
umn density of C3N - (1.6 xl0 12 cm -2 ) mat c hes th at deter- 
mined for IRC+ 1 0216 by IThaddeus et al.l i|200Sl ) exactly. 
ICernicharo et all l|2008h derived a rather low column den- 
sity for C 5 N - of 3.4xl0 12 cm -2 (lower than for CN - ), but 
note that this figure is likely to be an underestimate. We 
calculate a column density of A r (CsN - ) = 1.3xl0 14 cm -2 . 

At lower metallicity, small anions are a factor of a few 
lower in abundance than at Galactic metallicity. This is sur- 
prising given that the corresponding neutrals are more abun- 
dant, and the number of free electrons in the Magellanic 
envelopes is larger. The larger anions CeH - and C7H - are 
predicted to be more abundant at LMC metallicity than 
Galactic. 



4.3 Silicon chemistry 

Silicon chemistry was not included in the model of 
ICordiner fc Millarl (|2009l ). but is included here due to the 
potentially important reaction: 



C2H2 + Si 



SiC 2 



H 2 



(2) 



which could be a major sink of C2H2 in the inner regions of 
the envelope at low metallicity. However, the model shows 
that reactions with other hydrocarbons are the major loss 
mechanisms for C 2 H2, even in the highest metallicity case 
(i.e., where the abundance of silicon atoms is highest). 

SiC 2 was recently de tected by the Herschel Space Ob- 
servatory in IRC+10216 (|Cernicharo et al.ll2010l ). where its 
distribution traces the dust-for mation zone of the st ar, sim- 
ilar to SiO and SiS emission (|Fonfria et al.l 120081 ). How- 
ever, it also has a signifi c ant abundance in the outer 
envel ope (|Lucas et al.l Il995l ; iGensheimer, Likkel. fc Snvderl 
Il995h . as we find in our modelling (see Fig. |SJ), with 
the distribution peaking at ~5xl0 16 cm and reaching 
x(SiC2)=5x 10 -7 . These figures correspond well with the ob- 
servations of Lucas et al.l (|l995f ) and the chemical model of 
ICernicharo et al" 1 2010h . The abundance of SiC2 increases 
through the reaction between acetylene and silicon atoms 
(Eq. 0, and it becomes destroyed outside _R=10 17 cm by 
reaction with abundant C + ions. The peak abundance of 
SiC2 is similar in all metallicity regimes, and column den- 
sities are slightly enhanced in the LMC and SMC models 
(1.9xl0 15 :4.2xl0 15 :2.6xl0 15 ; MW:LMC:SMC). The chem- 
istry is slightly different in the LMC and SMC models, since 
the neutralisation of SiC2H + to form SiC2 dominates over 
reaction (J2]). SiC2H + itself is partially formed from SiC2, via 
the chain: 



Si + + C 2 H 



SiCf + H 
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Table 5. Column densities and ion-neutral ratios for anionic species in IRC+10216. 



Species 


Observcd column 


Neutral-ion 


Reference 




Calculated column 


Calculated 




density (cm -2 ) 


ratio 






density (cm -2 ) 


neutral-ion ratio 


CN" 


8.0X10 12 


400 


Agundez et al. f 20101 




3.0X10 12 


~8 000 


C 3 N- 


1.6xl0 12 


200 


Thaddeus et al. (20081 




1.6xl0 12 


~240 


C 5 N~ 


3.4xl0 12 


2 


Cernicharo et al. ("20081 




1.3xl0 14 


6 


C 2 H" 


<7xl0 10 


>70 000 


Agundez et al. (20101 




8.1xl0 8 


~10 8 


C 4 H" 


7.1X10 11 


4 200 


Cernicharo et al. (20071 




7.2xl0 14 


17 


C 6 H" 


4.0xl0 12 


16-100 


McCarthy ct al. (2006); Cernicharo et al. 


(20071 


8.2xl0 14 


4 


CgH" 


2.4xl0 12 


3-4 


Rcmiian et al. (2007); Kawaguchi et al. ( 


'20071 


2.1xl0 14 


5 



Table 6. Main silicon repositories at differing metallicitics. 



Species 


MW 


LMC 


SMC 


SiS 


70% 


65% 


64% 


SiO 


17% 


21% 


23% 


Si 


6% 


6% 


5% 


Si+ 


6% 


1% 


6% 


Siti 4 


2% 


2% 


2% 


SiC 2 


<1% 


2% 


2% 



C + + SiC 2 — ¥ SiCJ + c 
SiCl + H 2 — >• SiC 2 H+ + H 
SiC 2 H + + e" (or anion) — > SiC 2 + H ( + neutral) 

The major repository for silicon in circumstellar en- 
velopes is SiS (Table [6)l, in line with the high initial abun- 
dances of this molecule (Table [TJ. This is true in all metal- 
licity regimes tested, although the abundance of SiS drops 
slightly with metallicity as sulphur is preferentially incorpo- 
rated into CS. 



4.4 The effect of assumptions about physical 
parameters 

In order to fully understand the chemistry at low metal- 
licity distinct from the effect of the differing physical con- 
ditions, we calculated three models adopting the physical 
conditions of the Galactic model (Table [3]), and using the 
initial chemical abundances of the three different metallic- 
ity regimes. Hence, any variation in abundance or column 
density is solely due to the chemistry of the circumstellar 
envelope, or the initial abundances adopted. For clarity, we 
will name these the "fixed physics" models, as opposed to 
the "full physics" models, as described previously. 

In Table[7]we present the column densities derived from 
the full physics models, and their ratios. In the final column 
we have listed the ratio of column densities derived from the 
fixed physics models. Looking at how the two final columns 
in the table change gives us an indication of how the adopted 
physical conditions affect the chemistry of the model. 

For many of the species metallicity is the dominating 
factor, with column densities falling with metallicity despite 
an increase in the overall H 2 column density. Many of the 
nitrogen-bearing species are included in this group, which 
makes them good tracers of metallicity: CN, CH3CN, C5N - . 
However, there are exceptions to this: C3N, C5N, C3N - and 
the cyanopolyynes increase in relative column density at 



LMC metallicity, but decrease at SMC metallicity. Species 
containing heavy metals such as silicon or magnesium also 
reflect changes in metallicity, with the exception to this be- 
ing SiC 2 . For SiC 2 the effect of the lowering metallicity is to 
produce a decrease in the column density (see final column 
of Table [7)l . However, column densities in the full physics 
models increase for LMC and SMC metallicities, counter- 
acting the effects due to chemistry alone. Another species 
with similar properties is CS, which is produced inherently 
less in lower metallicity environments, but receives a boost 
in production in the denser environments of the full physics 
models. 

CO is not a very good tracer of metallicity in carbon 
stars, despite it being a good tracer of oxygen abundance un- 
der the fixed physics conditions. Observationally, CO emis- 
sion is often optically thick, and thus the observed flux in- 
tensity is not representative of the abundance of CO. In 
addition, in our models column densities of CO in the full 
LMC and SMC models are comparable, despite a factor of 
two difference in metallicity. CO is a very good tracer of the 
molecular envelope of carbon stars, though, and gives us an 
indication (see Fig. [5]) that in general, the envelopes of Mag- 
ellanic carbon stars will only extend ~70% of the distance 
of similar Galactic carbon stars. 

The column densities of hydrocarbons such as C 2 H 2 , 
C4H, C4H 2 and CeH 2 increase significantly at lower metal- 
licity, as we have seen. Table [7] shows that much of this in- 
crease is due to the changes in chemistry at lower metallicity, 
with the relative proportion of carbon increasing over oxy- 
gen (Table [2} . A significant increase in column density for 
these species is also due to the differing physical conditions 
found in Magellanic carbon stars. These species are formed 
through neutral-neutral reactions, and the higher densities 
of Magellanic envelopes means that these reactions proceed 
more quickly. This trait is not observed to the same degree 
in anions of hydrocarbon species, which are generally more 
abundant than in Galactic carbon stars for larger species 
(CeH~ and larger), but only by a factor of a few. This is 
somewhat surprising, since the ionisation fraction of the en- 
velope almost doubles from Galactic model to LMC model 
to SMC model. Smaller hydrocarbon species (C4H - , C5H - ) 
are more abundant in Galactic CSEs. 



5 OBSERVABLES 

With the impending completion of ALMA in 2013, it is inter- 
esting to consider the potential for ALMA to observe molec- 
ular line emission from extragalactic AGB CSEs, focusing 
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Table 7. Calculated column densities (cm -2 ) for species of interest. 



Species 


Nmw 


Wlmc 


NSMC 






Ratio 




Ratio fl f od . 

physics 


H> 


1 ssx 1 o 22 


5.77x 10 22 


1.15x 10 23 


I 


00 


2.00 


4.00 


I 


00 


1.00 


1.00 


C 2 H 


7 ssxin 16 


7 33xl0 16 


2.64xl0 16 


I 


00 


0.93 


0.33 


I 


00 


1.71 


1.78 


C2H 


8.12x10 s 


7.48x10 s 


l.lOx 10 9 


I 


00 


0.92 


1.36 


1 


00 


0.94 


1.02 


C'2 H2 


5.28xl0 17 


1.06xl0 19 


2.89x 10 19 


I 


00 


20.08 


54.77 


1 


00 


7.95 


8.76 


C2S 


2.67x 10 14 


3.43xl0 14 


1.78xl0 14 


I 


00 


1.28 


0.67 




00 


0.99 


0.50 


CH3CN 


1.65xl0 13 


1.88xl0 12 


2.81xl0 12 


I 


00 


0.11 


0.17 


1 


00 


0.18 


0.07 


CN 


2.46xl0 16 


1.48xl0 16 


8.84x 10 15 


I 


00 


0.60 


0.36 


I 


00 


0.62 


0.35 


CN _ 


3.01xl0 12 


2.93xl0 12 


1.46x 10 12 


I 


00 


0.97 


0.48 


I 


00 


0.92 


0.47 


CO 


1.29xl0 19 


9 93xl0 18 


1.04x 10 19 


I 


00 


0.77 


0.81 


1 


00 


0.39 


0.20 


cs 


5 66x10 16 


7.67xl0 16 


8.01xl0 16 


I 


00 


1.36 


1.42 




00 


0.65 


0.32 


C3H 


5.92xl0 14 


1.04xl0 15 


6.29x 10 14 


I 


00 


1.75 


1.06 


1 


00 


2.01 


1.81 




s nsxin 14 


7.37xl0 14 


8.19xl0 14 


I 


00 


2.43 


2.70 


1 


00 


2.80 


2.53 


C-aN 


3.92xl0 14 


6 OSxlO 14 


1.87xl0 14 


I 


00 


1.55 


0.48 


1 


00 


1.63 


0.95 


C 3 N~ 


1.64xl0 12 


1.77xl0 12 


7 93X10 11 


I 


00 


1.08 


0.48 




00 


1.13 


0.60 


C3S 


1.03xl0 15 


1.47xl0 15 


5.87xl0 14 


1 


00 


1.43 


0.57 


1 


00 


1.20 


0.62 


C4H 


1.21xl0 16 


2.00xl0 16 


fi 57x1 fl 15 


1 


00 


1.65 


0.54 


1 


00 


3.30 


3.86 


C4H~ 


7.18xl0 14 


3.04xl0 14 


l.llxlO 14 


1 


00 


0.42 


0.15 


1 


00 


0.99 


0.93 


v^4±±2 


2.89xl0 16 


4.35xl0 17 


3 55x1t) 17 




00 


15.04 


12.26 




00 


11.63 


12.85 


C5N 


7.51 xlO 14 


1.84xl0 15 


6.28xl0 14 




00 


2.45 


0.84 


1 


00 


2.79 


1.80 


CkN~ 


1.33xl0 14 


4.10xl0 13 


9.11x 10 12 


1 


00 


0.31 


0.07 


1 


00 


0.52 


0.29 


C 6 H~ 


8.24xl0 14 


1.84xl0 15 


1.73x 10 15 


1 


00 


2.23 


2.09 


1 


00 


2.16 


2.31 




1.15xl0 16 








00 


33.55 


30.50 


I 


00 


24.63 


29.99 


C 8 H" 


2.14X10 14 


7.35X10 14 


8.22X10 14 


1 


00 


3.44 


3.85 


1 


00 


2.07 


2.13 


HCN 


5.45X10 17 


1.44X10 18 


1.94X10 18 


1 


00 


2.64 


3.55 


1 


00 


1.20 


0.71 


HC 3 N 


6.44X10 15 


1.27X10 16 


5.26X10 15 


1 


00 


1.97 


0.82 


1 


00 


1.74 


1.02 


HC 5 N 


1.89X10 16 


8.66X10 16 


3.97X10 16 


1 


00 


4.58 


2.09 


1 


00 


3.59 


2.31 


HNC 


2.95X10 14 


9.11X10 14 


7.46X10 14 


1 


00 


3.08 


2.53 


1 


00 


4.77 


3.62 


MgNC 


3.92X10 15 


3.01X10 15 


l.SlxlO 15 


1 


00 


0.77 


0.46 


1 


00 


0.39 


0.15 


SiC 


l.lOxlO 14 


5.22X10 13 


3.02X10 13 


1 


00 


0.48 


0.27 


1 


00 


0.43 


0.24 


SiC 2 


1.93X10 15 


4.17X10 15 


2.56X10 15 


1 


00 


2.16 


1.33 


1 


00 


0.92 


0.47 


SiN 


1.05X10 14 


2.48X10 13 


1.19X10 13 


1 


00 


0.24 


0.11 


1 


00 


0.15 


0.04 


SiO 


4.37X10 16 


3.49X10 16 


3.70X10 16 


1 


00 


0.80 


0.85 


1 


00 


0.40 


0.21 


SiS 


1.77X10 17 


1.08X10 17 


1.03X10 17 


1 


00 


0.61 


0.58 


1 


00 


0.31 


0.15 



List of molecules observed bv lHe et al.l j2008l) in IRC+10216, with the addition of those in bold face. 



on the Magellanic Clouds. We have estimated line intensi- 
ties for nearly half a million lines assuming local thermal 
equilibrium (LTE) (Sect. [STTJ, and selected a few potential 
candidate lines for further investigation using a more rigor- 
ous non-LTE radiative transfer code (Sect. [572)) . 

5.1 LTE estimates 

Following the calculations of lOlofssonl |2008ah one can 
write: 



M 


1.2 


"15" 


1.6 




0.7 


' 1 " 


10~ 6 _ 








io- 3 _ 




D 



Jy, (3) 



where the CO (J=2-l) line flux density is given in terms 
of the mass-loss rate, expansion velocity, fraction of CO 
with respect to H2, and the distance, D, to the star. Insert- 
ing typical values for the LMC gives us an estimated flux 
density of ~0.04 Jy. Similarly, for the SMC, we calculate 
Sco(2-i) ~0.09 Jy, since we expect the density of circum- 
stellar envelopes in the SMC to be higher. Comparable flux 
densities are expected for the CO (J— 3-2) line. Such stars 
should be easily detectable within an hour's observations for 
the full ALMA array (5<r = 6 mjy at 2 km s _1 resolution). 
Stars with properties similar to our LMC stellar character- 
istics should be detectable out to «150 kpc within one hour, 



whilst stars similar to our SMC stellar characteristics could 
poten tially be seen ou t to distances of ~200kpc within an 
hour (Olofsson 2008a). 

For speci es other than CO , the line flux density can also 
be estimated (|Olofssonll20"08al ). For LMC carbon stars, 

e -£i/Mx 

S » 1.62 x 10- 9 gu A ul /xfi e n/rr - Jy, (4) 

using our standard parameters, and where g u , A u \ and E\ 
are the quantum mechanical degeneracy of the upper en- 
ergy level of the transition, the Einstein coefficient for the 
transition, and the lower energy level, respectively. Q(Tx) is 
the partition function, dependent on the excitation temper- 
ature of the molecule, Tx, w hich we assume to be 10 K for 
all rotational transitions (cf . IWoods et al.l 120051) . R c is the 
radius of the emitting region for each molecule, taken from 
the molecular distributions in the chemical model. Similarly 
for the SMC, 

e -BJkT K 

S « 3.72 x W- 9 g u A ul f x R e ——3y. (5) 

Estimated flux densities for molecular lines stronger than 
0.1 mjy are given in Table [5] The peak line intensities cal- 
culated from Eqs. [4] and [5] generally underestimate line 
strengths by a factor of ~5-50 in comparison to those from 
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Table 8. Predicted line strengths for carbon stars in the LMC, using Eg. l4l 



Band 3 (89-119 GHz) 



Band 6 (211-275 GHz) 



Band 7 (275-370 GHz) 



Molecule 


Frequency 
(GHz) 


Peak Flux 
(mjy) 


Molecule 


Frequency 
(GHz) 


Peak Flux 
(mjy) 


Molecule 


Frequency 
(GHz) 


Peak Flux 
(mjy) 


v-j4rl2 


89.315 


0.1 


SiO 


217.105 


0.7 


CS 


293.912 


0.3 


C4H2 


89.687 


0.1 


13 CO 


220.399 


1.0 


SiO 


303.927 


0.3 


SiC2 


93.064 


0.1 


CN 


226.632 


0.2 


13 CO 


330.588 


1.7 


C4H2 


97.834 


0.1 


CN 


226.660 


0.6 


CN 


340.008 


0.1 


CS 


97.981 


0.1 


CN 


226.664 


0.2 


CN 


340.020 


0.1 


C4H2 


98.245 


0.1 


CN 


226.679 


0.2 


CN 


340.032 


1.1 


C4H2 


98.655 


0.1 


CN 


226.874 


0.6 


CN 


340.035 


0. 1 


C4H2 


107.175 


0.1 


CN 


226.875 


0.9 


CN 


340.035 


0.7 


13 co 


110.201 


0.1 


CN 


226.876 


0.4 


CN 


340.248 


1.0 


CN 


113.491 


0.1 


CN 


226.887 


0.1 


CN 


340.248 


1.1 


CO 


115.271 


2.5 


CN 


226.892 


0.1 


CN 


340.249 


0.8 


SiC 2 


115.382 


0.2 


CO 


230.538 


Ri40 


CS 


342.883 


0.1 


C 4 H 2 


116.105 


0.1 


CS 


244.936 


0.5 


CO 


345.796 


^40 








SiO 


260.518 


0.5 


SiO 


347.331 


0.1 








C 2 H 


262.004 


1.5 


C 2 H 


349.338 


1.3 








C 2 H 


262.006 


1.1 


C 2 H 


349.339 


1.0 








C 2 H 


262.065 


1.1 


C 2 H 


349.399 


1.0 








C 2 H 


262.067 


0.7 


C 2 H 


349.401 


0.8 








C-C3H2 


265.759 


0.1 


HCN 


354.505 


58.9 








HCN 


265.886 


66.7 


C-C3H2 


368.294 


0.4 



Band 9 (602-720 GHz): Strongest lines include C2H3", HCN and CO, all with intensities on the order of lmjy. 



the non-LTE model ( £|5. 2|) . The exception to this is lines of 
HCN, which are over-estimated by ~9-16 times in compar- 
ison to the NLTE estimates. Our assumption of Tx = 10 K 
for HCN lines is probably a significant underestimate for 
a species which is only abundant in the inner regions of 
the circumstellar envelope. Adopting Tx = 75 K reduces 
Shcn(3-2) ^27 mjy, approximately halving the estimate for 
T x = 10 K. 



5.2 Non-LTE estimates 

In addition to our simple estimates for the expected line flux 
density described above, we have used the one-dimensional 
non-LTE Monte Carlo radiative transfer code, RATRAN 
l|Hogerheiide fc van der Takll2OO0l ). to estimate the line pro- 
files and line strengths of molecular line emission from 
evolved carbon stars in the LMC and SMC. In our simple 
estimates (see Eqs. [4] and [5}, we assume LTE and an excita- 
tion temperature of 10 K for all rotational transitions. In re- 
ality, the envelopes of AGB stars have a radi ally-dependent 
temp erature and density structure (see e.g.. lJura fc Morrisl 
Il98ll ) which limits the accuracy of calculations made with 
these assumptions. The temperature varies from ~1000K 
in the inner envelope to ~10K in the outer envelope. Simi- 
larly, the density structure exhibits a R~ 2 behaviour so that 
in the outer envelope the density is much lower than the 
critical density of rotational transitions (A u i/T,iC U i, where 
i indicates all energy levels lower than the upper level, u). 
We used molecular data from the Leiden Atomic and Molec- 
ular Database (LAMDA0) which has tabulated energy lev- 
els, Einstein A coefficients and collisional rates for many 



molecules with transitions in the (sub-)mm region of the 
electromagnetic spectrum l|Schoier et al.ll2005l ). 

In our radiative transfer calculations, as in the physi- 
cal and chemical modelling, we treat the stellar envelope as 
spherically symmetric with a constant outward stellar wind 
velocity of 10 km s" 1 for the LMC and 5 km s" 1 for the SMC. 
Since our line profiles are broadened to a width >20kms 
and >10 kms -1 , respectively, we use a spectral resolution of 
lkms -1 for both sources. We assume a distance to source 
of 50kpc for the LMC and 66kpc for the SMC. We gen- 
erated line profiles and line strengths for rotational transi- 
tions of CO, CS, CN, HCN, SiO and SiS which fall into the 
expected spectral range of ALMA "Full Science" operations 
(~30 GHz to «950 GHz). These are the molecules which our 
calculations suggested may be sufficiently abundant to pos- 
sess emission strong enough to be observable and for which 
collisional rates are available. Using the ALMA sensitivity 
calculatojf], we determined those line transitions which may 
be observable with ALMA Full Science within a realistic 
observing time. Our results are listed in Table [9] We de- 
termine that molecular transitions in Bands 6 and 7 make 
particularly good targets. Assuming a spectral resolution of 
2 kms -1 and an array size of 50 antennae, we find the obser- 
vation times for LMC and SMC molecular transitions range 
from a few minutes for the CO J=2-l and J=3-2 transi- 
tions at 230 and 345 GHz to several hours for transitions 
of HCN, CS, SiO, SiS and CN which also fall into Bands 
6 and 7. In addition, the CO J=l-0 transition, in Band 3, 
may also be observable to high signal-to-noise within one 
hour in evolved stars within both galaxies. We note here 
that our time calculations are likely upper estimates since 
the ALMA full array may eventually consist of up to 66 



http : //www. strw. leidenuniv.nl/~moldata/ 



almascience . eso . org/call-f or-proposals/sensivity- calculator 
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Figure 7. Synthetic spectra of an LMC carbon star for ALMA Band 6 and Band 7. 
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Figure 8. Synthetic spectra of an SMC carbon star for ALMA Band 6 and Band 7. 



antennae in total. Also, the online ALMA sensitivity cal- 
culator overestimates the predicted observing times relative 
to the more accurate calculator available using the ALMA 
Observing Tool, again, indicating that our predicted times 
are likely upper estimates. In Figs. [7]-[8] we present the ex- 
pected line spectra due to emission from the listed molecules 
from a carbon star in both the LMC and SMC, for ALMA 
bands 6 and 7. Line profiles of each line transition are also 
provided in online-only material (Figs. 9-32). Reducing the 
mass-loss rate of the models by a factor of three reduces the 
line strengths in Table [5] by a similar factor. Thus probing 
the molecular inventory of carbon stars with low mass-loss 
rates in the MCs will be challenging. 

Looking firstly at the predicted line spectra for the LMC 
(Fig. [7)), we see that the strongest transitions are due to CO, 
CN, HCN and SiO. Multiple transitions of CN, SiO, CS and 
SiS are available within a single band which would allow de- 
termination of the temperature and density of the emitting 
gas, enabling some constraints on the physical conditions in 
the envelope. For the SMC (Fig. [8}, we see a similar pattern, 
however, we now find that transitions of SiS are comparable 
in strength with those of SiO and HCN. Comparing the two 
metallicity regimes, we see differences in the line strength 
ratios between CO/CN and SiO/SiS, with these ratios gen- 



erally decreasing with metallicity. This follows, given the de- 
creased amount of oxygen available in the SMC relative to 
the LMC. We also see the line strength ratios of CN/HCN 
slightly increase with metallicity, likely related to the in- 
creased strength of the ISRF in the SMC. Our calculations 
demonstrate that there are observable tracers of metallicity 
and physical conditions in extragalactic carbon stars which 
ALMA will allow us to probe. 



6 DISCUSSION 

This investigation was prompted by the observation of 
very common and very deep molecular absorptions due to 



C2H2 in the mi d-infrared spectra ( 5-38 /jm) of LMC car- 
bon stars (e.g., IWoods et al.l l201ll . for a number of ex- 
amples) compared to Galactic carbon stars. This occurs 
particularly in the extreme carbon stars, which are los- 
ing mass at very high rates. Of course, the fact that 
these deep features ex i st ha s been established for sev- 
eral years: ISpeck et al.l (|2006h analys ed the deepes t C2H2 
absorption feature observed to date; IZiilstra et al.l (|2006l ) 
noted that acetylene b ands were stronger at low metal- 
licity; iLeisenring et al. J2008I) measure d eq uivalent widths 
of molecular lines; Sloan et alj {2006) and lLagadec et al.l 
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LMC SMC 



Molecule 


Transition 


Frequency 
(GHz) 


Peak Flux 
(mjy) 


S/N 


Integration Time 


Peak Flux 
(mjy) 


S/N 


Integration Time 


CO 


J 


= 1-0 


115.271 


19.2 


10 


50.4 min 


19.8 


10 


47.2 min 




J 


= 2-1 


230.528 


64.2 


10 


1.3 min 


65.3 


10 


1.3 min 




J 


= 3-2 


345.796 


69.0 


10 


2.8 min 


67.5 


10 


3.0 min 




J 


= 4-3 


461.041 


65.4 


10 


2.6 hr 


57.8 


10 


3.9 hr 




J 


= 6-5 


691.473 


64.2 


5 


1.4 hr 


52.9 


5 


2.5 hr 


HCN 


J 


= 3-2 


265.886 


4.1 


5 


1.6 hr 


2.6 


5 


4.2 hr 




J 


= 4-3 


354.505 


6.7 


5 


1.7 hr 


1.2 


5 


4.5 hr 


CS 


•J 


= 5-4 


244.936 


2.8 


5 


2.6 hr 


2.7 


5 


2.6 hr 




J 


= 6-5 


293.912 


3.8 


5 


2.1 hr 


3.5 


5 


2.8 hr 




.J 


= 7-6 


342.883 


4.7 


5 


2.7 hr 


1.3 


5 


2.8 hr 


SiO 


J 


= 5-4 


217.105 


3.9 


5 


1.4 hr 


2.7 


5 


3.6 hr 




J 


= 6-5 


260.518 


5.0 


5 


1.0 hr 


3.6 


5 


2.0 hr 




J 


= 7-6 


303.927 


5.9 


5 


1.0 hr 


4.3 


5 


1.8 hr 




J 


= 8-7 


347.331 


6.4 


5 


1.3 hr 


1.8 


5 


2.4 hr 


SiS 


.1 = 


= 12-11 


217.817 


2.9 


5 


2.5 hr 


2.7 


5 


3.6 hr 




.1 = 


= 13-12 


235.961 


3.1 


5 


2.7 hr 


3.0 


5 


2.8 hr 




.1 = 


= 14-13 


254.103 


3.3 


5 


2.0 hr 


3.3 


5 


2.0 hr 




.1 = 


= 15-14 


272.243 


3.4 


5 


2.1 hr 


3.5 


5 


2.2 hr 




.1 = 


= 16-15 


290.380 


3.5 


5 


2.6 hr 


3.6 


5 


2.7 hr 




.1 = 


= 17-16 


308.516 


3.5 


5 


3.2 hr 


3.7 


5 


3.3 hr 


CN 


J 


= 2-1 


226.333 


2.3 


5 


3.6 hr 


1.9 


5 


5.7 hr 




J 


= 2-1 


226.659 


6.0 


5 


37.6 min 


1.1 


5 


1.1 hr 




J 


= 3-2 


340.031 


4.4 


5 


2.7 hr 


3.7 


5 


4.6 hr 




J 


= 3-2 


340.249 


9.1 


5 


39.9 min 


6.9 


5 


1.1 hr 



In our integration time estimates, we assume an array size of 50 antennae and a spectral resolution of 2.0km/s. 



detected acetylene absorptions in SMC carbon stars 
that appeared deeper than in corresponding Galactic and 
LMC carbon stars. There have also been studies in the 
acetylene-d ominated 3.1/3 . 8 fim b a nds which show simi- 
lar re s ults: van Loon et all (Il999al). Matsuura et all (|2002l . 
120051 ). Ivan Loonl l|2006j ). Ivan Loon et all l|200Sl ) performed 
a 3 /»m band study, finding that acetylene absorption fea- 
tures in LMC and SMC carbon stars were equivalent in 
depth. The origin of the 3/im acetylene b and is a mixtur e 
of photospheric and circumstellar origin (|van Loonl l2006h . 
whereas the 13.7 fim absorptio n seen in Spitzer sp ectra 
is predominantly circumstellar (|Matsuura et al.ll2006t ). We 
have shown via thermal equilibrium calculations that "pho- 
tospheric" acetylene increases in abundance as the metal- 
licity is lowered (Fig. [1]), qualitatively matching what we 
observe at 3 fim. We have also shown via radially-dependent 
non-equilibrium calculations that indeed acetylene is more 
abundant in carbon-rich circumstellar envelopes at lower 
metallicity (Fig. [2]| . Initial investigations into carbon stars 
in other low(er)-metallicity gala xies show that ace tylene fea- 
tures are also very strong (e.g.. Isioan et al1l2009l ) Unfortu- 
nately we do not yet have sufficiently good quality infrared 
data to perform quantitative analyses. 

Related to the study of acetylene is the study of hy- 
drogen cyanide. It absorbs in the 3.1 fj,m band with acety- 
lene, and also forms part of the 13.7pm band, absorbing at 
~14/xm. There is also a band at 7/im, close to the acety- 
lene band at 7.5 /im. Its contribution to these bands has 
been detected in a number of Magellanic c arbon stars, but 
not nearly so prevalently as acetylene (e.g . , iMatsuura et al.l 
120051 ; IZiilstra et al"1l2006j ; Ivan Loonl 12009 ). In some objects 



the absorptions are strong (|Matsuura et al. 20021). whereas 
in oth er objects, absorptions are very weak ( Matsuura et al.l 
2008). Both the thermal equilibrium model and the circum- 
stellar envelope model predict high abundances of HCN at 
low metallicity, despite the lower abundance of elemental ni- 
trogen. This would imply that both the 3 /im and 13.7 /jm 
bands should show ample evidence of HCN absorption. Our 
high abundances may be an effect of our choice of 3M© 
nucleosynthesis models (Sect. 12. ip since the nitrogen abun- 
dance in AGB stars is dependent on (initial) stellar mass, 
such that stars with high lumi nosity (i.e., high i nitial mass) 
should have abundant HCN (|Matsuura et all [2005). The 
variation in feature strength seen in the infrared may re- 
flect varying initial stellar mass in the observed samples. 
The impact of different choices of stellar mass should be 
looked into in future models. Scaling initial abundances up 
or down uniformly generally results in a similar scaling of 
CSE abundances; however, changes in stellar mass are un- 
likely to produce such a uniform scaling, but would result in 
abundance enhancements in some elements over others. 

Emission from two large circumstellar molecules was ob- 
served recently in Magellanic objects. The fullerenes Ceo 
and C70 ha ve been detected in a decade of planetary neb- 
ulae (PNe; ICarcfa-Hernandez et al.l l201ll ; IZhang fc Kwokl 
l201ll provides a summary of all fullerene detections) 
in the LMC, as well as several Galactic objects, in- 
cluding (proto-)planetary nebulae and reflection nebulae. 
The Magellanic clouds are presumably rich in fullerenes, 
whether they form from smaller hydrocarbons coagulating 
or th e destruction of hydrogenated a morphous carbon dust 
(e.g., iGarcia-Hernandez et alj l201ll ). In Galactic nebulae 
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fullerenes potentiall y take up one pe r cent of the elementa l 
carbon abundance (Cami et al ] |20ld ; ISellgren et a"l]|2010l ); 
as Table[2]shows, there is significantly more carbon available 
for fullerene formation in Magellanic objects. On average, 
iGarcfa-Hernandez et alj 's sample of Magellanic carbon-rich 
PNe contains 60% more Ceo than the Galactic PN, Tc-1. As 
a representative of large molecules, C23H2 is two orders of 
magnitude more abundant in Magellanic carbon stars than 
in Galactic (Fig. [3} in our models. As a representative of 
cyclic molecules, benzene is ~200 times more abundant in 
Magellanic carbon star envelopes than Galactic; its column 
density is iY(C 6 H 6 )~10- 5 Af(C2H 2 ). 

These four molecules currently make up the extent of 
our knowledge of gas-phase species in the envelopes of car- 
bon stars in the Magellanic Clouds. Given the derived line 
intensities and simulated spectra from Tableland Figs. [7]- 
[8] ALMA will be able to detect a handful of molecules in a 
reasonable amount of time: detections of CO rotational lines 
will enable us to derive mass-loss rates and envelope expan- 
sion velocities; detections of other molecules will enable us to 
verify our chemical models, and to improve our understand- 
ing of circumstellar carbon chemistry at low metallicity. 



7 SUMMARY 

Driven by the recent interest in deep molecular absorption 
features in Spitzer Space Telescope spectra of evolved car- 
bon stars in the Magellanic Clouds, we have conducted the 
first investigation into the circumstellar chemistry of carbon 
stars in three metallicity environments: that of the Galaxy 
(solar metallicity), the LMC (half solar metallicity) and the 
SMC (one fifth solar metallicity). The general trend is that 
the abundances of hydrocarbons are greatly enhanced at low 
metallicity, so much so that the larger members of hydrocar- 
bon families are more abundant than smaller members (e.g., 
C6H2 and C8H2 are more abundant than C4H2). At just half 
of solar metallicity, acetylene becomes more abundant in 
the circumstellar envelope than carbon monoxide, the most 
abundant molecule apart from H2 in Galactic CSEs. This is 
indicative of a suppressed oxygen chemistry at lower metal- 
licity: other oxygen-bearing species are also less abundant. 
Nitrogen chemistry is also suppressed, except when it is in- 
corporated into hydrocarbons (e.g., cyanopolyynes). This 
means that the main nitrogen repository in the envelope 
shifts from N2 in Galactic stars to HCN in Magellanic stars. 
A similar trend is seen for heavy metals (e.g., Si), where 
SiC2 behaves like a hydrocarbon, but other silicon-bearing 
molecules are rarer. 

Physically, the lower expansion velocities of Magellanic 
CSEs mean that the envelopes are more dense. This implies 
that neutral-neutral reactions in the inner envelope are more 
rapid, but the stronger radiation fields in the LMC and SMC 
cause the molecular envelopes to be lesser in extent than in 
Galactic CSEs. The ionisation fraction increases as metallic- 
ity drops, and the main charge carrier moves from a heavy 
metal (Mg + ) in Galactic envelopes to C + . 

We have compared the results from the Galactic model 
to the carbon star IRC+10216, and achieve reasonably good 
agreement in terms of molecular abundances and distribu- 
tions. We discuss the qualitative results that have been ob- 
tained from infrared observations of carbon stars in the Mag- 



ellanic Clouds, and also achieve good agreement, although 
constraints are few. Finally, we use our models to predict 
line intensities and images as seen by ALMA. ALMA is ex- 
pected to improve our understanding of carbon stars in the 
Magellanic Clouds, and lines of CO, HCN, CS, CN, SiO and 
SiS should be detectable within 3 hours of observation in 
LMC carbon stars. Lines of C 2 H and 13 CO are predicted to 
be somewhat weaker, but may also be observable. 
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